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Abstract. The behavior of the aromatic features as a function of metallic- 
ity/radiation field hardness in M101 HII regions and starburst galaxies is inves- 
tigated using Spitzer/IRS spectra. These spectra reveal that the 6.2, 7.8+8.6, 
and 11.3 am. aromatics have different dependencies on metallicity /radiation field 
hardness. Specifically, the 6.2 and 7.8+8.6 am. aromatics are weak or absent at 
a metallicity, 12+log(0/H), of ^8.0 while the 11.3 /im feature remains relatively 
strong. These results apply to both the M101 HII regions and starbursts show- 
ing that HII regions can indeed be thought of as "mini-starbursts." Comparison 
of this work with various candidate materials results in a good match with the 
annealing behavior in Quenched Carbonaceous Composite (QCC) which is an 
amorphous, hydrogenated, carbonaceous solid. 

The aromatic features are a family of dust emission features seen in many 
dusty astrophysical environments. The strongest aromatics are seen at 3.3, 
6.2, 7.7, 8.6, 11.3, 12.7, and 17.1 /xm and all the aromatics have been iden- 
tified with C-H and C- C bending and stretching modes of hydrocarbons con- 
taining aromatic rings ( Tielensl 2005f ) . The sh ape and strength of the aro- 
matics have been seen to vary in single sources (Werner et al.ll2004bh . 



among 



various Galactic so urces J van Diedermoven et all hoo4 and "arn^ galaxies 
( Engelbracht et al.l 1200,4 iMadden et al.ll2006l : IWu et al.ll200fil L A number of 
different materials have be en proposed as th e carriers of these features since 



they were first discov ered dGillett et al 



phous Carbon (HAC. iDulev h Williams! _ , v 

posites (QCC. ISakata et al.l ll984\ Polycyclic Arom atic Hydrocarbons (PAHs 



1973) including Hydrogenated Amor- 
1983), Quenched Carbonaceous Com- 



Allamandola et alll985T ). coal ((Papoular et alll989h , and nanodiamonds ( Jones d'Hendecourtl 
2000). The leading candidate material for the aromatics is PAH molecules and 
the y are included in the modified "astrono mical PAHs" form in dust grain mod- 
els (|Li fc Draineil200ll : IZubko et al.l E004 \. The preception that PAHs are the 
carrier of the aromatics can be seen by the common use of the term "PAH 
features" to refer to the aromatics. The identification of PAH molecules with 
the aromatics may be somewhat premature given that laboratory spectrum of 
a mix of PAHs doe s not match the observed wave lengths and strengths of the 
aromatics in detail (Hudgin s Sz Allamandolall2004h . 

As part of the Multiband Imaging Photometer for Spitzer (MIPS) Guar- 
anteed Time Observations (GTO), two programs were carried out to explore 
the variations in the infrared (IR) Spectral Energy Distributions (SEDs) of 
massive star formation as a function of metallicity. The first program targets 
starburst galaxies with metallicities [12+log(Q/ H)1 between 7.2 to 8.9 , with 
emphasis on the lowest metallicity starbursts (|Engelbracht et aDl2005| K The 
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Figure 1. The HII regions with IRS spectra are identified on the MIPS 
24 jum image. The HII region Hodge 681 is also identified (but not observed 
with IRS) as one of the lowest metallicity HII region detected in the MIPS 
images. 



second program targets M101 as it is a large , face-on spiral galaxy with HII re- 
gions which have metallicities from 7.4 to 8.8 (Zari tskv et al Has 4; v an Zee et al.l 
1998: lKennicutt et al.ll2003l ) . The combination of these two programs probe HII 
regions, the components of starbursts, and entire starbursts. 

Understanding the IR SEDs of starbursts and HII regions is important to 
understanding dust properties as well as the SEDs of redshifted galaxies under- 
going massive star formation. Probing how dust emission features vary with 
physical environment (e.g., metallicity) can give clues to the dust grain mate- 
rials. Interpreting number counts in deep fields to probe the star formation 
history of the Universe requires a good understanding of galaxy SEDs; the be- 
havior of the aromatics directly affects models of number counts. For example, 
the aromatics moves through the MIPS 2 4 um. band at redsh ifts between 1 and 
2, causing a peak in the number counts ( Laeache et al.ll2~004l h 



Observations 

The starburst and M101 HII regions have bee n observed with all the instru- 
ments on the Spitzer Space Telescope (Spitzer IWerner et al. I l2004al h This pa- 
per concentrates o n preliminary results from the Infrared Spectrograph (IRS, 
iHouck et al.ll2004bl). Prelimin ary results f or M101 from the I nfrared Array Cam- 
er a (IRAC. iFazio et aTll2004l ) and MIPS (|Rieke et al.ll2004l ) images were given 
bv lGordon et al.H2005lh Preli minary results from the IR AC and MIPS imaging 
of the starbursts were given bv lEntrelbracht et al.l ()200.^ . 
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Figure 2. IRS spectra of M101 HII regions are shown normalized at 10 /jm 
in the top plot and the same spectra are shown in the bottom plot offset by 
an arbitrary amount. The full IRS spectra extend from 5-38 /im. The HII 
region names are listed along with th eir measured metallicities and excitation 
temperatures l|Kennicutt et alJl2003|) . 



The spectra of the M101 HII regions (Fig. ^) were taken in the spectral 
mapping mode with full slit width steps perpendicular to the slit. Because 
the ShortLow and LongLow observations were taken approximately 3 months 
apart, the orientations of the slits were approximately the same. These spectral 
mapping observations were reduced with the SSC pipeline (version S12.0.2) and 
combined using CUBISM (Smith, J.D.T. 2006, in prep.). Spectra were extracted 
from the cubes for each HII region using an aperture with a radius of 15" and 
sky annuli from 20 to 30" referenced at 24 //m. The size of the aperture and sky 
annuli were varied linearly with wavelength to roughly account for the changing 
diffraction limited PSF. The resulting spectra are shown in Fig. [21 
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Figure 3. IRS spectra of starburst galaxies are shown normalized at 10 /im 
in the top plot and the same spectra are shown in the bottom plot offset by 
an arbitrary amount. The full IRS spectra extend from 5-38 /im. The galaxy 
names are listed along with their measured metallicities (Moustakas, J., in 
prep.). 



The spectra of the starburst galaxies were taken in staring mode. The 
reductions were done using the SSC pipeline (version Sll.0.2). The redundant 
observations were combined and the spectra were extracted with a combination 
of SPICE and IRAF. The extractions were for the full slit, since the galaxies 
were usually noticeably extended. The final spectra have been re-sampled to a 
common wavelength scale resulting in some degree of smoothing. The resulting 
spectra are shown in Fig. |31 
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Figure 4. The equivalent widths of 3 aromatics are plotted on the left for 
the M101 HII regions. The ratio of the aromatic to the 11.3 /J.m aromatic is 
shown on the right. The measurements for the M101 nucleus are plotted with 
open symbols. 



Results 

The M101 HII region and starburst spectra shown in Figs. El & El clearly display 
the behavior of decreasing strength of the aromatics with metallicity which is 
easiest to see in the strong 7.8+8.6 /im aro matic complex . Thi s is expected 
from previous Infrared Space Observatory (iMadden et al.l 12006) and Spitzer 
( Houck et aTll2004al : lEneelbracht et al.ll2005l : IWu et al.ll2006T ) observations. The 
spectra presented here allow us to determine how individual features vary as a 
function of metallicity. From visual inspection of the spectra, the 6.2, 7.8, and 
8.6 /iin features are weak or absent around a metallicity of 8.0 as seen from the 
strength of these features in NGC 5471 (Fig. EJ) and UM 462 and Tol2138-405 
(Fig. |31). This is not the case for the 11.3 and 12.7 /im features which are seen in 
all the spectra in both plots except for the lowest metallicity starburst spectra 
which are too noisy. Thus, different aromatic features seem to display different 
behaviors with metallicity. 

We have quantified the behavior of the aromatics in the M101 HII region 
spectra by simply measuring the equivalent widths and total flux in the 6.2, 
7.8+8.6, and 11.3 /um aromatic features/complexes. The measurements were 
done for these three features/complexes as it is possible to identify continuum 
on either side of them. In addition to these aromatics, we also measured the 
fluxes in the [S IV] 10.6 /im and [S III] 18.7 /im lines as the ratio of these two 
atomic emission lines is a measure of the hardness of the radiation field. 
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The equivalent widths and flux ratios of these aromatics are plotted in Fig.Q] 
normalized to Searle 5, the strongest aromatic emitter, versus the [S IV] /[S III] 
ratio and metallicity. The weakening of the aromatics seems better correlated 
with increasing radiation field hardness ([S IV] /[S III]) than metallicity. The 
plots of the aromatics to the 11.3 /mi aromatic show that the 11.3 fj,m aromatic 
weakens at a much slower rate than the 6.2 and 7.8+8.6 /im aromatics. In 
addition, the changing of the aromatic ratios also seems better correlated with 
radiation field hardness than metallicity. 

This preliminary analysis suggests that the weakening of the aromatics is 
due to the hardening of the radiation field and that the aromatic spectrum 
changes as the features weaken. The pattern of aromatic spectru m change is 
not i n agreement with the predictions for PAH molecule ionization ( Bakes et alJ 
120011 . Fig. 5) which retain the 6.2+7.8+8.6 fim complex at higher ionizations, 
but predicts weak or no 11.3 /xm aromatic. The aromatic spectrum change 
does seem to be in qualitative agreement with the annealing behavior of QCC 
( Sakata et al~lll984l . Fig. 2). 

A more complete analysis of the M101 HII regions will be published by Gor- 
don et al. (in prep.) including a more accurate decomposition of the aromatics 
using the PAHFIT program (Smith, J.D.T., in prep.). A more complete analysis 
of the starburst spectra will be published by Engelbracht et al. (in prep.). 
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